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Abstract — This paper is concerned with the problem of design- 
ing computationally efficient Generalized Comb Filters (GCF). 
Basically, GCF filters are anti-aliasing filters that guarantee 
superior performance in terms of selectivity and quantization 
noise rejection compared to classical comb filters, when used as 
decimation filters in multistage architectures. 

Upon employing a partial polyphase (PP) architecture pro- 
posed in a companion paper, we develop a sensitivity analysis in 
order to investigate the effects of the coefficients' quantization 
on the frequency response of the designed filters. 

We show that the sensitivity of the filter response to errors 
in the coefficients is dependent on the particular split of the 
decimation factor between the two sub-filters constituting the PP 
architecture. The sensitivity analysis is then used for developing 
a fixed-point implementation of a sample filter from the class of 
GCF filters, used as reference filter throughout the paper. 

Finally, we present computer simulations in order to evaluate 
the performance of the designed fixed-point filters. 

Index Terms — CIC-filters, comb, decimation, decimation filter, 
delta, delta-sigma, fixed-point, GCF, generalized comb filter, 
partial polyphase, polyphase, EA, sigma, sigma-delta, sine filters. 



I. Introduction and problem formulation 

The design of computationally efficient decimation filters 
for oversampled SA A/D converters [l]-[3], as well as for 
classical oversampled A/D converters, has received a renewed 
interest recently, spurred by intense research activities in con- 
nection to the design of digital front-ends for both wideband 
digital receivers and Software Defined Radio receivers [4] -[6]. 

Consider a base-band analog signal x(t) (with bandwidth 
+f x }) sampled by an Analog-to-Digital (A/D) converter 
at rate f s = 1/T„ = 2pf x > 2f x , where p(> 1) is 
the so-called oversampling ratio. If p is close to unity, the 
A/D converter operates at the Nyquist frequency, whereas 
for p » 1 we are referring to oversampled A/D converter. 
When p 1, the decimation of the oversampled discrete-time 
signal x(n/ f s ) is usually accomplished by cascading two (or 
more) decimation stages, followed by a FIR filter that provides 
the required selectivity on the discrete-time signal x(nTj^) at 
baseband. 

Fig.Q]shows a multistage decimation architecture composed 
by m decimation stages operating on the oversampled signal 
x(nT s ). Also shown are the data rates of the sampled data at 
the input, as well as at the output, of each decimation stage 
in the multistage decimation architecture. 

For the sake to contain the computational complexity of the 
overall architecture, the first decimation stage usually employs 
a multiplier-less filter [7]-[8]. A widely used filter featuring 
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this property is the comb filter [2], [3], [9], which provides 
an intrinsic anti-aliasing effect by placing its zeros in the 
middle of each folding band, i.e., in the integer multiples of 
the digital frequency 1/D (D is the decimation factor). The 
transfer function of a A c -th-order comb filter is defined as [2]: 

/II — • V 

H c (z) = 



1 1 



D 1 



(i) 

where D is the decimation factor. 

Unless the design of classical FIR filters, the design of 
a decimation filter embedded in a multistage architecture 
imposes stringent constraints on proper bandwidths. Let us 
elaborate. 

Consider the architecture shown in Fig. [1] where the over- 
sampling ratio p is factorized as 



P : 



Up, 



In the previous relation, each Pi is a proper positive integer. 
The sampling rates at the input and output of the i-th stage 
are, respectively, 



fi-i = fi ■ Pi, Vi = 1, . 



with f 



fs, and 

fi = 



Vi = 1, 



, m. 



With this setup, consider the frequency response Hi(fd) of 
the i-th decimation filter pictorially shown in Fig. [2] The 
digital frequency is the normalized signal bandwidth at 
the input of the i-th decimation filter. Notice that, for i — 1, 
it is f° = f x / f s = (2/o) _1 : this is the normalized bandwidth 
of the signal sampled by the A/D converter at rate f s . For any 
other i, the relation between /* and f° is 

f i c = ft 1 P i , Vi = l,...,m. 

Given this setup, the frequency response Hi(fd) has to atten- 
uate the quantization noise (QN) within the frequency bands 
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because the QN falling inside these frequency bands, will fold 
down to baseband (i.e., within the useful signal bandwidth 
[—fc^ 1 , +/c _1 ]) due to the sampling rate reduction by Pj in 
the i-th decimation stage [10]. Such a QN will irremediably 
affect the signal resolution after the multistage decimation 
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Fig. 1, General architecture of a multi-stage decimation architecture for oversampled A/D converters. Given p, the discrete-time signal x(n/f B ) is decimated 
through a cascade of m decimation stages, obtaining the Nyquist-sampled discrete-time signal ic(nTjv) at baseband. The i-th decimation stage is composed 
by an anti-aliasing filter Hi(z) followed by a decimator by Pi. The filter Hi(z) operates at the data rate fi—\, while the i-th decimation stage delivers a 
discrete-time signal with rate /j = /j_i/P£, The index i takes on values in the range {!,..., m}. 



architecture. On the other hand, the frequency ranges labeled 
as don't care bands in Fig. |2j do not require stringent selec- 
tivity, since the QN within these bands will be rejected by the 
subsequent anti-aliasing filters in the multistage chain. 

With this background, let us provide a survey of the recent 
literature related to the problem addressed in this paper. 
Tutorials on the design of multirate filters can be found in [11], 
[12], while essential books on this topic are [7] -[8]. The design 
of optimized multistage decimation and interpolation filters 
has been recently addressed by Coffey in [13]-[14], while 
the design of multistage decimation architectures relying on 
constituent cyclotomic polynomial filters has been presented 
in [10]. A 3-rd-order modified decimation sine filter was 
proposed in [15], and developed in [16]. The class of comb 
filters was generalized in [17], whereby the author proposed an 
optimization framework for deriving the optimal zero rotations 
of GCFs for any filter order and decimation factor D. 

Other works somewhat related to the topic addressed in 
this paper are [18]-[24]. In [18] and [19], the authors pro- 
posed computational efficient decimation filter architectures 
for implementing non recursive classical comb filters. In [20], 
the authors proposed the use of decimation sharpened filters 
embedding comb filters, whereas in [21]-[22] the authors 
addressed the design of a novel two-stage sharpened comb 
decimator. In [23], the authors proposed novel decimation 
schemes for SA A/D converters based on Kaiser and Ham- 
ming sharpened filters, then generalized in [24] for higher 
order decimation filters. Papers [25]-[27] focus on the design 
of decimation filters with improved frequency responses. 

The main aim of this paper is to develop a mathematical 
framework for the design of fixed-point GCF decimation filters 
relying on a partial polyphase FIR architecture proposed in the 
companion paper [28]. To this end, we develop a sensitivity 
analysis in order to investigate the effects of the coefficients' 
quantization on the frequency response of the designed filters. 

For conciseness, we focus on the design of the first decima- 
tion filter in the multistage architecture in Fig. Q] even though 
the proposed analysis can be easily extended to the design of 
the other anti-aliasing filters in the cascade. 

The sensitivity analysis paves the way to a statistical ap- 
proach useful to identify the coefficient word lengths of the 
proposed fixed-point architecture. Moreover, we show that the 
proposed analysis highlights some key issues in connection 
to the choice of the proper split of the decimation factor 
between the polyphase stage and the cascaded FIR sections 
of the employed partial polyphase architecture. 

The rest of the paper is organized as follows. In Section ITfl 
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Fig. 2. Pictorial representation of the frequency response of the i-th 
decimation filter, Hj(z), along with the key frequency intervals to be carefully 
considered during the design. P$ is the decimation factor, is the data 

rate at the input, /j is the data rate of the decimated data, and fl~ 1 is the 
normalized signal bandwidth. 

we briefly review the transfer functions of GCF filters, as well 
as the partial polyphase architecture employed throughout the 
paper, and outline the key advantages that these filters feature 
with respect to classical comb filters. Section [III] presents a 
mathematical framework for evaluating the sensitivity of the 
frequency response of GCF filters to the quantization of the 
coefficients. In Section ITVl we discuss general guidelines for 
the design of the proposed filters, and present some simulation 
results. Finally, Section [VI draws the conclusions. 

II. Overview of GCF filters: The Partial 
Polyphase Architecture 

Our objective in this section is to recall the fundamental 
concepts for analyzing both the time-domain behaviour and 
the frequency response of GCF filters, as well as to highlight 
the main differences between GCF and classical comb filters. 

For conciseness, our discussion will be restricted to a 3rd- 
order GCF filter, which will be used as a reference scheme 
throughout the paper, and we will present the non-recursive, 
partial polyphase architectural! developed in the companion 
paper [28]. 

'Even though both recursive and non recursive implementations can be 
devised for GCF filters, the non-recursive architecture does not present any 
instability problem deriving from the quantization of the coefficients. We 
notice in passing that recursive GCF filter realizations rely on zero-pole 
cancellations, which can be impaired by the quantization of the coefficients. 
This is the basic reason for the use of this FIR architecture in the developments 
that follow. 
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Fig. 3. Zero locations of GCF niters within the k-th folding band. Zeros are displayed in the z-plane, as well as in the frequency domain, in order to 
highlight the main differences between GCF and classical comb filters. 



Let us focus on the design of the i-th decimation filter, Hi(z), 
in Fig. Q] and, for ease of notation, assume P,; = D and omit 
the subscript i in Hi(z). Given D, and recalling the definition 
of the folding bands in (O, a classical 3rd-order comb filter 
(see ([TJ with 7V C = 3) presents 3-rd-order zeros in the complex 
locations 

z k = e^ k , Vfc = l,...,D-l, 

or, equivalently, in the digital frequencies f Zk = jy, k € 
{1, . . . , fejvf }. Therefore, a 3rd-order zero is placed in the 
middle of each folding band. This idea is illustrated in Fig. [3] 
where the /c-th folding band is shown in the z-plane, as well 
as in the frequency domain. 

On the other hand, a 3rd-order GCF filter places, in the £>th 
folding band, 3 zeros in the digital frequencies ~ — and 



z-transfer function of the non recursive filter decimating by 
£>2- The latter function is defined as 



-2' 



Hn{z) = Yl [ 1 + r i- \ z - +z 

i=P P + l 

whereby the coefficients r{ are defined as 

n = 1 + 2 cos (2*a) , Vi = p p + 1, 



2-2 1 



+ Z 



-3-2 1 



,P-1, 



(4) 



(5) 



whereas a = 2 • 0.79 • irf^ 1 . We notice in passing that 
according to the commutative property in [9], the filter Hjsi{z) 
can be realized with the cascade of p — p p — 1 stages, each 
one decimating by 2. The integer p p can take on any value in 
the set {-1, . . . ,p— 1}. 

After some algebra, the frequency response of the filter 



£r, whereas the edges of the fc-th folding band are 4 — Hn(z) can be evaluated by substituting z 



D 

p~ x and +/c . Therefore, as shown in Fig. [3] the choice 
a = q2nfl , with q 6 [0, +1], allows a better distribution of 
the three zeros within the k-th folding band, whose width is 
strictly related to the bandwidth / c i_1 of the useful discrete- 
time signal. 

The optimal parameter q has been found in [17], and we 
will employ such a value throughout this work. As an example, 
the optimal value q = 0.79 is such that a 3rd-order GCF filter 
features an additional 8dB of QN rejection over a classical 
3rd-order comb filter around the folding bands. 

Once again, let us focus our attention on the 3rd-order 
GCF filter, and consider a decimation factor D that can be 
expressed as the p-th power-of-two, i.e., D — 2 P , where p is a 
suitable integer greater than zero. Moreover, let us factorize the 
decimation factor D as D = D\ ■ D2, whereby D\ = 2 Pp+1 , 
£>2 = 2 p ~ Pp ~ 1 , and p p can take on any integer value in the 
set {— l,...,f> — 1}. With this setup, the z-transfer function 
of a third-order GCF filter realized with the partial polyphase 
FIR architecture in Fig. |U can be defined as follows: 



H N {e* u ) 



-J3-2 1 



(6) 



(7) 



H(z) = H ■ H P (z) ■ H N (z), 



(3) 



whereby H Q is a constant term ensuring unity gain at base- 
bancQ. The function Hp(z) is the z-transfer function of the 
polyphase section decimating by D\, whereas Hjy(z) is the 

2 For simplicity, we omit this constant term in the derivations that follow. 
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• [cos (3 • 2 i - 1 w) + n cos (2 l " 1 w)] , 

whereby uo — lirfd, and r, is defined in (O. 

The impulse response hp(n),yn € [0, 3Z?i — 3], whose z 
transfer function is identified by Hp(z), is defined as [28]: 

h P (n) = e +]an e~ 2jaH J e JQfc2 ^ x t( k i 

fc 3 =0 k 2 =0 fei=0 

The definition of the sequence xt (n) in (Q is 

x t (n) = 5(n)-r5(n-Di)+rS(n-2D 1 )-5(n-3D 1 ), (8) 

whereby r = 1 + cos(aZ?i) and a = q2-Kp ( T 1 . 

The impulse response in (O is used to obtain the polyphase 
components 

e*(n) = h P (Dm + k), Vfc e [0,D a - 1] (9) 

of the filters Ek(z) in the architecture shown in Fig. |4] 

Let us spend a few words about the parameters noticed in 
H(z). The choice p p = p — 1 allows the GCF filter H(z) to 
be fully realized in polyphase form, whereas the value p p = 
— 1 is such that the filter is realized as the cascade of p non 
recursive decimation stages, each one decimating by 2. Any 
intermediate value of p p S {0, ... ,p — 2} yields the partial 
polyphase decomposition depicted in Fig. |4] 
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iH,(z) 



D, E (z) 



D, 



E,(z) 
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Fig. 4. Architecture of the partial polyphase implementation of the decimation filter H(z). 



The first polyphase decimation stage allows the reduction of 
the sampling rate by D\, thus reducing the operating rate of 
the subsequent decimation stages belonging to Hn(z). Any 
stage of Hn(z) in Fig. |4]is constituted by a simple FIR filter 
operating at a different data rate. Such an example, the i-th 
stage, with i £ {0, . . . ,p — p p — 2}, is characterized by the 
transfer function \l + ri (z^ 1 + z~ 2 ) + z~ 3 ] operating at rate 



follows: 



converter. 



where f s is the sampling frequency of the A/D 



III. Design of Fixed-Point GCF Filters 

In this section we consider the problem of evaluating 
the sensitivity of the filter H(z) in (01 with respect to the 
coefficients enclosed in both Hp(z) and Hn(z). Then, the 
sensitivity function is employed in a design algorithm that 
defines statistically the size of the registers in the fixed-point 
implementation of the filter in such a way that the error 
function in the frequency domain between the filter iJ(e J ") 
and the filter H q (e? u ), which employs quantized coefficients, 
is within given bounds with a preassigned probability. 

We show that the proposed framework gives a precise 
answer on the choice of the proper split of the decimation 
factor D between the two substages Hp{z) and Hm{z). It 
is anticipated that the best solution from a sensitivity point 
of view consists in implementing the filter H(z) without the 
polyphase stage Hp(z), i.e., with D = D% and D\ = 1. 

Before proceeding further, let us derive some observations 
on the sensitivity function employed throughout this section. 
Given a frequency response H (e JW ), the sensitivity analysis 
is usually accomplished on the magnitude of H (e- JW ) with 
respect to its coefficients. However, the sensitivity analysis 
based on the use of the frequency response H (e 3 ") can be 
derived much easier than the one that employs the function 
^(e^)! [29]- Moreover, the sensitivity function related 
to H (e JW ) provides an upper bound to the one related to 
\H (e 3 ") . Consider the frequency response 



Hi 



and a tagged multiplier m belonging to H < 



Then, the 



dm 



a M ei ")\ j>p(H(e>")) 



dm 

-j \H(ei" 



dm 



Upon observing that both \H (e juJ ) \ and (p (H (e juJ )) are real 
functions of w, the bound 



dH (e 



dm 



> 




\H\ 



V dm 



(10) 



easily follows. Owing to this result, the sensitivity function 
used in this work relies on the derivatives of the frequency 
response H (e JW ) with respect to its multipliers. 

When the multipliers belonging to H (e^) are quantized by 
employing rounding, the magnitude of the frequency response 
becomes: 



\H q {e juJ ) \ = \H(e>")\+A\H(e ju )\ 



(ID 



derivative of H (e^j with respect to m can be evaluated as 



whereby A Lff(e JW )| is an error function that measures the 
distortion of the ideal frequency response iJ(e ja; )| from the 
one employing quantized coefficients. Recalling the definition 
of the folding bands given in (O, a key observation in the 
proposed framework is that the error function A|ff(e ja; )| 
must be properly bounded only within the folding bands, 
concisely identified by FB. Therefore, care must be devoted 
to the sensitivity analysis only within the folding bands, 
while the behaviour of the error function outside the FB does 
not affect the proposed fixed-point design. Notice that these 
considerations only hold for the design of decimation filters in 
multistage architectures, and cannot be extended to the design 
of classical FIR filter. 

The aforementioned considerations can be formalized as 
follow: 

| A \H{e iu )\\ < x(w),Vw e FB, (12) 

whereby is a suitable-positively defined- tolerance func- 
tion. Even though we can theoretically choose to differentiate 
the behaviour of the tolerance function x(w) among the 
various folding bands, we do not pursue this approach in this 
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work. Therefore, the functions x( w ) use d i n the following will 
be constant functions across the folding bands. 

Next line of pursuit consists in investigating a statistical 
technique in order to identify the word-lengths of the filter 
coefficients. To this end, assume that the frequency response 
H(e : ' ul ) contains N coefficients rounded by employing the 
same fixed-point resolution (6-bit rounding quantization). Af- 
ter quantization, each multiplier rtii can be written as 

m g< i = m l + Ami. 

Therefore, the quantization of the N coefficients yields 
N zero-mean, statistically independent random variables 
Ami that are identically and uniformly distributed in 
[-2- 2b /2,+2-' 2b /2] [33]. Under these hypotheses, and re- 
membering the relation (fTTT i. the variance of the error function 
A (^(e 3 ")) can be evaluated as follows: 



°A|ff(e^)| 



2 

°Ar 



N 

,£ 

i=l 



dH(e juJ ) 



dm; 



2 

a Ar, 



Sr(e 



3U\ 



(13) 



whereby o\ m = 2~ 2b /12 is the variance of the random 
variable Am under the hypothesis to employ rounding to the 
nearest quantization level. 

Owing to the condition JV > 1, the error function 
A Lff(e 7W ) can be modeled as a zero-mean Gaussian random 
variable with variance given by ( TT31 [31]-[32]. Therefore, we 
can estimate the probability 



P 



P[\A\H\\<ya A]H{ 



e*")|J 



that A | (e JtxJ ) | falls within a proper interval, say from 

-y<TA\H(ei«)\ t0 +2/0-A|H(eJ")|. as follows: 



-yCTA|ff( e j") 



2na 



A I H< ^'I dx. 



(14) 



It is convenient to employ the new variable 

x 



V2, 



CA|-ff(e3"")| 



,2 



= -= e~ z dz. (15) 



in ( TPfl i. thus obtaining 

p = P[|A|F|| <ya AmeJ u 

Let us spend few words about the result (15[ . The term p 
is the probability that the magnitude of the error function 
| A l-H^e-^)] | in (fT2l is upper-bounded by y<JA\H{e^)\- 

The relation between p and y in (TT~5T > is illustrated in Fig. [5] 
As an instance, the value y — 2 has to be chosen in order 
to guarantee with a probability equal to 95% that the error 
function is bounded by 2a A \H{ei")\ m tne frequency domain. 

How can we employ this result in a practical design? Upon 
recalling (fT2l . and given a proper p, we choose y in such a 
way that the following relation holds: 



y ■ 0-A|iT(ei«)| » x(u) => y ■ O-A, 

This relation can be rewritten as 

2' ' 



(16) 



1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

y 

Fig. 5. Relation between the probability p and the parameter y in fT5}. 

By doing so, we guarantee that (fT2l is verified with probability 
p given by (fT~5T >. 

Given this statistical framework, let us focus on the fixed- 
point design of the considered GCF filter. 

Assume that the filter coefficients are represented with the 
following fixed-point notation: /„ bits are devoted to the 
integer part of the coefficients, while F n is the number of 
bits devoted to the fractional part. Therefore, the size of the 
filter coefficients is equal to 1 + /„ + F n bits, accounting for 
the sign of the number. 

The next two subsections derive the sizes of both I n and 
F n in the fixed-point implementation. 

A. Evaluation of the fractional size, F n 

Considering b = F n , and solving (fTot for F n , we can obtain 
the size F n of the fractional part in order for (fT2l to hold with 
a probability p given by ST5t : 



F„ 



log 2 



12 min 



,gFB yy/S T {ei") 



(17) 



whereby the minimum is taken only over the folding bands 
derived in (0, and ["•] is the ceil of the underlined number. 

The evaluation of the fractional part F n in ( TTTb relies on 
the sensitivity function St ( e "'")- To keep the presentation 
concise, the derivation of the sensitivity function is reported 
in the Appendix. 

The behaviour of F n in (TTTb as a function of the decimation 
factor D\ of the polyphase stage is illustrated in Fig. [6] for 
various values of the function x( w )- The setup for deriving 
the results in Fig. [6] is as follows. We considered a two-stage 
decimating architecture (to = 2 in Fig.[TJi and an oversampling 
ratio p = 4 • D, where D is the other parameter associated to 
each subplot in Fig. [6] The first decimation filter in the 2-stage 
architecture is the investigated GCF filter, thus Hi{z) |,=i = 
H(z) in Fig. [2 

The normalized digital bandwidth of the useful signal at the 
input of the GCF decimation filter is f^ 1 | »=x = f° = j^- 
This is also half the width of the folding bands seen by the 
first decimation filter in the two-stage architecture. Therefore, 
the value of a appearing in the definition of filters Hp(z) and 
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Fig. 6. Behaviour of F n in (T7J as a function of the decimation factor D\ of the polyphase stage for three values of the function x(u)), assumed constant 
across the folding bands. Left subplots are associated to the value y = 2 in \\5\ corresponding to the probability p = 95%, whereas the rightmost subplots 
are associated to the value y = 1.63, which corresponds to the probability p = 90%. The three curves in each subplot are parameterized with respect to the 
value of the function x(uj) as follows: curves labeled with the mark — X — are for x(a>) = 5 X 10 -3 , the mark — * — is associated to x( £i) ) = 10 — 3 , and 
mark — * — is used for x(a>) = 10 — 4 . 



H N (z) in © and Q is 0.79 • 2tt/° = 0.79w/p. As a note 
aside, notice that the value of a would be different if the GCF 
filter were used in the second stage of the multistage chain in 
Fig. Q] due to the different value of 

We considered three different constant functions^ x( w ) m 
order to draw (jTTJ, namely x( w ) = 5 x 10~ 3 , x( w ) = 10~ 3 , 
and xH = 1(T 4 

Each subplot is associated to a specific decimation factor 
D = Di ■ D 2 - Therefore, given the constant D noticed in 
the ordinate of each subplot, the number of stages belonging 
to Hm{z) is reduced as long as D\ increases. In particular, 
the abscissa D\ = 1 in each subplot is associated to the case 
D = L>2, which means that the GCF filter is implemented 
without the polyphase stage. 

Moreover, we consider two different values of y related to 
the probabilities p — 90% and p = 95% illustrated in Fig. [5] 

Some observations are in order. 

• A comparison among the leftmost and the rightmost 

3 We notice in passing that other behaviours can be associated to the 
tolerance function x( w )> depending on the required sensitivity desired in 
the various folding bands. However, we do not pursue this approach in this 
work, and assume that all the folding bands affect equally the size F n of the 
fractional part. 



subplots in Fig.|6]reveals the need of one additional bit for 
the fractional part in order to guarantee that the constraint 
on the tolerance function is attained in the frequency 
domain with probability 95% with respect to the case 
90%. 

• For given y, D, and D\ in the abscissa, the number F n 
of fractional bits increases as long as a lower tolerance 
function x( w ) is desired. 
> Given D, the size F n of the fractional part of the fixed- 
point implementation increases as long as the number of 
cascaded cells in H^{z) decreases. This is equivalent to 
say that F n increases as long as Di does. This obser- 
vation suggests that the GCF filter H(z) implemented 
as Hn(z), i.e., without the polyphase stage, allows to 
contain the computational complexity of the GCF filter. 
The latter observation above suggests that an effective imple- 
mentation of the GCF filter is H(z) = Hn(z). Therefore, the 
3rd-order GCF filter is realized with the cascaded architecture 
shown in Fig. [7] 

This architecture follows from (0]i upon setting p p = — 1: 



p— * 

Hn(z) = {J [l + r/t- (V 2 * + z- 2 - 2 *) 



+ z 



-3-2* 



18) 



k=t) 



7 



k-th decimation stage 




\+r p . l {z}+z 2 ) +z~ 3 




Fig. 7. Architecture of a non-recursive implementation of the decimation filter H(z). Index k, which identifies the decimation stage, belongs to the range 
0, 1, . . . , log 2 D — 1 = p — 1. The figure also shows an effective implementation of the fc-th filter with transfer function 1 + r^z -1 + z~ 2 ) + z~ 3 . 



whereby the coefficients r k are defined as 

r k = 1 + 2 cos (2 fc a) , Vfc = 0, . . . ,p - 1 
a = 2 • 0.79 • 7T f^ 1 . (19) 

Applying the commutative property [9], the cascaded imple- 
mentation shown in Fig. [7j easily follows. 

B. Evaluation of the integer part, I n 

This section is focused on the evaluation of the size of the 
integer part /„ in the fixed-point representation of the filter 
coefficients. To this goal, consider the architecture in Fig. [7j 
and focus on the fc-th decimation stage. Let 1% be the size of 
the integer part in the fc-th decimation stage. 

The impulse response associated to the transfer function 
1 + r k (z~ l + z~ 2 ) + z~ 3 is 

h k (n) = 6{n) + r k S(n - 1) + r k S(n - 2) + 5{n - 3). 

Upon relying on general considerations about dynamic range 
overflow, it is simple to observe that the worst-case dynamic 
range growth G k of the fc-th stage is 

G k < log 2 \h k (n)\) = log 2 (2 + 2r fc ) < 3, 
Vn=0 / 

where the last inequality stems from the observation 

r k = l + 2cos(2 fc a) < 3, Vfc = 0, . . . ,p - 1. 

Therefore, the size of the integer part 2* (in bits) that avoids 
overflow, is equal to the sum between the width of the input 
word (in bits) and G k - 

IV. Simulation Results 

In this section, we compare some GCF filters designed 
with the framework proposed in the previous section, with 
classical comb filters. We also provide a set of simulation 
results obtained by employing a fixed-point realization of a 
3rd-order GCF filter for decimating a discrete-time signal 
oversampled by a EA A/D converter. 



The first set of results is proposed to compare the frequency 
response of a 3rd-order GCF filter H(e^) with the one 
obtained with a fixed-point implementation. Let us summarize 
the setup. We consider a two-stage decimating architecture 
(to = 2 in Fig. [TJ, whereby the first stage employs a GCF 
filter decimating by D = D2 — 16 (i.e., the GCF filter is 
implemented with the architecture shown in Fig.|7jwith p = 4), 
whereas the second stage presents a decimation factor equal 
to 4. With this setup, the oversampling ratio is p = 64. 

From the upper-leftmost subplot in Fig. [6] we notice that 
F n — 7 in order to satisfy the bound 

\A\H{e^)\ \ < xH = 1(T 4 ,Vcj e FB, 

with a probability equal to 95%. 

The normalized digital bandwidth of the useful signal at the 
input of the GCF decimation filter is 

ri—l I ro ^ ^ 

/c |i=1 ~ Ic ~ 2p ~ 128' 

The magnitude of the frequency response (dotted-line curve) of 
the 3rd-order filter _ff(e- J ") without coefficients' quantization is 
shown in Fig.[8]for D = 16. From (O, we notice the presence 
of the following folding bands: 

k 1 . _k_ 1 1 1 l. -1 u C 

16 128 ' 16 "r 128 J ' ~ ' 1 ' ' ' M ~ °> 

some of which have been highlighted in Fig. [8] In the 
same figure, we show for comparison the frequency response 
(continuous curve) of the GCF filter whereby the coefficients 
have been quantized with F n = 7, as discussed above. Notice 
that the two frequency responses are mostly superimposed, 
thus confirming the effectiveness of the proposed design 
framework. 

The frequency response H q (e^) is compared with the one 
of a classical 3rd-order comb filter in Fig. [9] The figure clearly 
highlights the behavior of the GCF filter across the folding 
bands: unless a classical 3rd-order comb that places a 3rd- 
order zero in the frequencies ^g, Vfc = 1, . . . , fcj\/ = 8, the 
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|H(e>")|-dB 

- |H y-)l-6B 




Folding bands 




Fig. 8. Magnitude of the frequency response of filter H{e ]ul ) employing 
real coefficients (continuous curve), and frequency response of filter H q (e 3u} ) 
employing fixed-point coefficients (dotted curve). Both frequency responses 
are mostly superimposed. The decimation factor is D = f6. 



Fig. 9. Magnitude of the frequency response of a classical 3-rd-order comb 
filter H^{e J ^) (dotted curve), and frequency response of filter H q (e ]u} ) 
employing fixed-point coefficients (continuous curve). The decimation factor 
is D = 16. 



GCF filter places its zeros in the frequencies 

0.79 -J"; Vfc = l, 



A. A - 79 - r - 

16' 16 /c ' 16 



k M = 



The last set of results is obtained by resorting to simulation. 
Employing Matlab, we simulated a 2nd order EA converter 
with a 2-level quantizer and a sampling frequency f s = 
25.6 kHz. The input signal is a band-limited random signal 
with bandwidth f x = 100 Hz. From the values of f s and f x , 
it is p = 128, while the normalized digital bandwidth of the 
sampled signal is f° = ^ = ggg. 

The oversampled signal is then decimated by D = 16 
employing a 3rd-order GCF filtefl The power spectrum of 
the digital signal at the output of the EA A/D converter is 
shown in the upper subplot of Fig.QI)] Notice that, as expected, 
the useful signal with bandwidth f° = ^ is shrunk at 
baseband, while the EA A/D converter has pushed the noise 
power spectrum outside the useful signal bandwidth [0,/°]. 
The power spectrum of the decimated signal is shown in 
the lower subplot of Fig. |T0j Notice that the useful signal 
bandwidth is now =/ c °-C = |k 0.063. 



V. Conclusions 

This paper focused on the design of computationally effi- 
cient Generalized Comb Filters (GCF), i.e., anti-aliasing filters 
that, employed as decimation filters in multistage architectures, 
guarantee superior performance in terms of selectivity and 
quantization noise rejection compared to classical comb filters. 
GCF filters can be realized by relying on both IIR and 
FIR architectures, even though FIR schemes do not present 
instability problems stemming from coefficients' quantization. 

As a reference filter in the class of GCF filters, a third order 
FIR architecture, realized by employing a partial polyphase 
architecture, was used throughout the paper. We proposed a 
sensitivity analysis in order to first investigate the effects of 

4 We recall that the order of the GCF filter has to be greater or equal to 
B + I, whereby B > 1 is the order of the SA modulator [24]. 




- S (f a ) [dB] 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 





0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Fig. 10. Power spectrum of the digital signal at the output of a 2nd-order 
SA A/D converter sampling at / s = 25.6 kHz (upper subplot), and power 
spectrum of the signal decimated by D = 16 with a 3rd-order GCF filter 
(lower subplot). 



the coefficients' quantization on the frequency response of the 
designed filters, and, then, to define the registers 'lengths in the 
proposed fixed-point implementation. 
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Appendix 

In this Appendix we derive the sensitivity function 
St (e^) needed in the evaluation of the fractional part F n 
in (V7\ . Let us evaluate the function St (e? u ) m EJl for 
the considered partial polyphase architecture. To this end, 
we consider three different cases depending on the value of p p . 

First Case: p p = — 1. This is the case in which D\ = 1 
and D = D2- The frequency response of the GCF filter is 
H (e J ") = Hn (e JW ), and the filter is implemented with- 
out the polyphase stage. Therefore, the sensitivity function 
St (& u ) can be evaluated as follows: 



S T (e j ») = \H N (e> 



p-i 

i=0 



cos (3 • 2 2 



cos(2 i - 1 w) 



-2 

(20) 



The previous equation stems from (0 upon noting that, for 
p.p = —1, the derivative of H^{e : > u ) with respect to r u , Vu = 
0, ... ,p — 1, can be evaluated as: 

dH N {e ju 



nr4 [cos (3. 2 

By multiplying and dividing for the function 

cos (3 • T^uS) + r u ■ cos (T^uS) , 
and recalling ©, OTT) can be rewritten as follows: 

d JM^l = H N {en 



2e -j3'2 "C0S(2"- 1 ^)- 



r m • cos (2 m - 1 Lu)] 
(21) 



cos( 2 U 



cos(3-2 u 1 uj)-\-r u -cos(2 u_ ' 



(22) 



Second Case: p p = p — 1. This is the case in which D2 = 1 
and D — D\. The frequency response of the GCF filter is 
H (e JW ) = Hp and the filter is fully implemented with 

a polyphase architecture. 

The polyphase decomposition of the z-transfer function 
Hp(z) is defined as follows: 



From (l23T > and (O, Hp(eJ UJ ) can be rewritten as: 

Dl -i Li/OiJ 

H P (e luJ ) = J2 hp ( Dl ' 71 + fc)e~^ (Dl " +fc) , 

fc=0 n=0 

which is valid Vn, k such that < D\ ■ n + k < L. 
Upon observing that 



(23) 



(24) 



dH P {eP u ) 



dhp 

the sensitivity St reduces to: 

2 



St (e^) 



L 

E 



dH P 



= L = 3D 1 - 2. 



dhp(i) 

which corresponds to the number of multipliers in Hp{e : ' U3 ). 
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Intermediate Case: p p = 0, ...,p — 2. This is the case 
where the GCF filter is implemented with the partial polyphase 
architecture discussed above. Let N = L + p — p p — 1 be the 
number of multipliers belonging to H (e- 7 ' w ) (L = 3Z?i — 2 
is the number of multipliers belonging to Hp (e- J "), while 
p— p p — 1 is the number of coefficients belonging to H N (e- 7 ")). 

The sensitivity function St assumes the following 

expression: 



N 

E 



dH (< 



9mi 



El* 



27V 

i=l 

P-Pp-1 



drrii 

E i^i 2 



(25) 



i=l 



A' 



9m; 



S T (e*") = L-|H JV (e J ' w ) | 



p-p P -i 



After some algebra, the previous relation can be rewritten as: 

cos(3-2P"+ l - 1 w) 

CI 

i=l 

L\H N {e j ")\ 2 + \H(e j ") 
cos (3 • 2Pp +l - 1 u;) 



\H N (e j n\ E 



cos(2Pp+ i - 1 uj) 

|2 



Pp+j 



(26) 



p-Pp-i 

E 

i=l 



whereby 



cos(2Pp+ J - 1 w) 
r Pp+l = 1 + 2 • cos (2^ +l a) . 



